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A new parameterisation describing the relationship between viscosity and solid fraction valid 

O 

' ^ at large solid content was recently proposed in Costa [2005]: 
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^ ■ where r\ is the relative viscosity, <p the volume fraction of particles, B the Einstein coefficient 
(with a theoretical value B = 2.5), and a, (3, 7 are three adjustable parameters. The proposed 

in ; 

parameterisation approximately reduces to the classical, well established relationships for small 

O ' 

fractions. On the other hand, a weak point of equation (1) is that for large 0, the relative viscosity 
tends quickly to a constant value, because the non-linear term in the erf function rapidly saturates 
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it as approaches unity. 

For very large solid fractions, no reliable experimental data are available and there are also some 
intrinsic problems even in properly defining relative viscosity [see Costa, 2005]. However, 
we can reasonably assume that the trend of the viscosity-crystal content relationship obtained 
from data by van der Molen and Paterson [1979] at high solid fractions for partially-melted 
granite is generally valid, and this it can be roughly described by a power law relationship. In 
fact, there is no reason to assume that effective relative viscosity of melted rocks tends towards 
a plateau region as — > 1 as predicted by (1). From a quantitative point of view it is no 
simple to determine correctly the controlling parameters, such as the critical fraction at which a 
rheological transition occurs or the maximum value that effective relative viscosity can reach as 
— > 1 . Moreover there are several indications that these values strongly depend on the particle 
shape and particle size distribution. However, since magma viscosity controls magma transport, 
when modelling volcanic processes such as dome growth there is a strong need to estimate the 
viscosity dependence on crystal-content, even in the limit of very high crystal content where no 
experimental observations are available [see e.g., Melnik and Sparks, 2005]. 
Here we propose a simple modification of relationship (1) which improves all the main positive 
features of it, yet at the same time, does not show a plateau region as — > 1, generally increases 
the performance of the model, and exactly recovers the classical well established relationships 
valid for small fractions. In order to obtain a satisfactory description of the effective relative 
viscosity on the entire range of 0, from zero to near the unity, one more parameter is necessary, 



DRAFT 



February 2, 2008, 6:44am 



DRAFT 



COSTA: CRYSTAL-CONTENT VISCOSITY DEPENDENCE X - 3 

i.e. four parameters. The final four parameter model we adopt is the following: 




where < a < 1,0* represents the critical transition fraction, 7 > 1 is a measure of the rapidity 
of the rheological transition, and 8 controls the increase of 77 as </> — > I. We can note that, 
for large 5, the second term at the numerator represents a negligible correction when <f> < </>*, 
whilst it becomes important when (ft > </>*• It easy to see that as <f> < (/>* with decreasing </>, the 
parameterisation (2) tends exactly to the Krieger and Dougherty [1959] relationship: 




and as (ft — > 0, it recovers exactly the Einstein equation: 

77(0) Ml + B 4>) (4) 

where B is the Einstein coefficient. 

An example of the good performance of relationship (2) is shown in Figure 1 . Here we plotted 
relative viscosities deduced from the data of van derMolen and Paterson [1979] in the high crys- 
tal fraction regime, considering a crystal-free fluid viscosity of 10 4 Pa s as suggested by van der 
Molen and Paterson [1979]. These values are matched with and mean values of relative viscosi- 
ties reported by Thomas [1965]. Although we must keep in mind all the intrinsic limitations of 
these data [see Costa, 2005], we can observe that over the entire range of solid fraction, the model 
is able to reproduce all data with an excellent correlation (R 2 = 0.999) and the curve behaviour 
suggested by van der Molen and Paterson [ 1 979] . Another example is shown in Figure 2 where, 
as in Costa [2005], the relative viscosities for Mg3Al 2 Si 3 0i2 and Li 2 Si 2 5 measured by Lejeune 

andRichet [1995] and the fitted values are shown. However a direct comparison with Figure 2 
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of Costa [2005] is not possible because the relative viscosities values reported in Figure 2 of 
Costa [2005] were erroneously underestimated in systematic way for both Mg 3 Al 2 Si30i2 and 
Li 2 Si 2 5 . 
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Figure 1. Effective relative viscosities deduced from data by van der Molen and Paterson 
[1979] at high solid fraction (circles), and from mean value of relative viscosities reported by 
Thomas [1965] at low solid fraction (squares). For data by van der Molen and Paterson [1979], 
for crystal-free fluid viscosity at 800°C and 300 MPa, we considered a value of = 10 4 Pas. 
The correlation coefficient is 0.999. The best fit parameters are a = 0.999916, = 0.673 and 
7 = 3.98937, 6 = 16.9386. 
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Figure 2. Relative viscosities from data by Lejeune and Richet [1995] (triangles) and values 
using eq. (2) (line) for Mg 3 Al2Si30i 2 (left) and Li 2 Si 2 5 (right). The best fit parameters are a = 
0.9703; 0.9995 , 0* = 0.503; 0.409 and 7 = 0.924; 4.214, 5 = 40.386; 1.149for Mg 3 Al 2 Si 3 0i 2 
and Li 2 Si 2 5 respectively. 
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